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Abstract. We attempt to model Quark Gluon Plasma (QGP) evolution from the initial 
Heavy Ion collision to the final hadronic gas state by combining the Glauber model initial state 
conditions with eccentricity fluctuations, pre-equilibrium flow, UVH2+1 viscous hydrodynamics 
with lattice QCD Equation of State (EoS), a modified Cooper- Frye freeze-out and the UrQMD 
hadronic cascade. We then evaluate the model parameters using a comprehensive analytical 
framework which together with the described model we call CHIMERA. Within our framework, 
the initial state parameters, such as the initial temperature (Tinit), presence or absence of 
initial flow, viscosity over entropy density (?7/s) and different Equations of State (EoS), are 
varied and then compared simultaneously to several experimental data observables: HBT radii, 
particle spectra and particle flow, /nd^ values from comparison to the experimental data for 
each set of initial parameters will then used to find the optimal description of the QGP with 
parameters that are difficult to obtain experimentally, but are crucial to understanding of the 
matter produced. 



1. Introduction 

Recreating Quark Gluon Plasma (QGP), matter that existed in the first millisecond after the Big 
Bang, has been the goal of experiments at RHIC and now at the LHC. QGP is created during the 
course of a collision of massive relativistic ions, such as gold ions in the Brookhaven Laboratory 
or lead ions at CERN, and has been shown to exhibit the properties of a viscous liquid [lj[2]. 
To understand the properties of the matter created in these collisions, the evolution of matter 
during the collision itself must be understood. Currently, we understand the basic stages of 
matter evolution, but lack the nuanced understanding of interdependencies and variation in 
parameters required for precise description of the collision. 

Our present understanding of a heavy ion collision can be summarized as follows. After 
collision time r ^ 1 fm/c, a strongly interacting liquid of quarks and gluons is formed, which 
is described with various accuracy by hydrodynamical models. As the gas cools, mesons and 
hadrons "freeze out" and start interacting as colorless objects, until a phase of a relatively stable 
gas of hadrons is reached. 

In this work we attempt to evaluate numerically the validity of models describing this picture 
by systematically varying model parameters and calculating distributions of proximity of 
fit to the experimental data. The evaluating algorithm used for these purposes is called 
CHIMERA, which stands for Comprehensive Hydrodynamical Integrated Modelling Evaluation 
and Reporting Algorithm, and will be described in detail in the next section. 
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Figure 1. Hydrodynamical initial conditions energy density distributions for Au+Au collisions 
at y/sNN — 200 GeV with b =4.4 fm, scaled as Npart (cl, b, c), or Ncoii (d, e, /). Panels a and d 
were produced using optical Glauber geometry, panels b and c using Monte Carlo Glauber, and 
panels c and / using Glauber geometry with event-by-event fluctuations. 



2. Model description 

CHIMERA is hydrodynamics-based framework, thus we first start by determining initial inputs 
and an initial energy distribution. In the analysis presented here, the initial conditions are 
computed using Monte Carlo Glauber model of the colhding nuclei [3] (e.g.. Figure [l] panels b 
and e). However, CHIMERA is capable of handling several initial condition implementations, 
including optical Glauber initial conditions (Figure [l] panels a and d), Glauber conditions 
with fluctuations (Figure [Tj panels c and /), and Colour Glass Condensate-inspired initial 
conditions (not shown). We implement initial conditions with either the number of collisions 
(Ncoii) or number of participants (Npart) scaling. In the future, we plan to introduce a dial-up 
mechanism for this scaling, so that the initial conditions would scale with a predetermined mix 
of Ncoii and Npart, as is closer to experimental observations. The Equation of State (EoS) used as 
a CHIMERA input is obtained from Lattice QCD calculations by the HotQCD collaboration [5j, 
and agrees well with s95p-vi EoS by Huovinen and Petreczky This EoS is much softer than 
the resonance gas EoS predicted by the MIT Bag Model 0, but is stiffer than the lattice-inspired 
EoS used in VH2 Hydro [8J. All four equations of state are shown in Figure [2j 

As the last step before starting viscous hydrodynamics, we employ pre-equilibrium flow, i.e., 
flow that originates from the non-equilibrated matter at early times (approximately before r = 1 
fm/c), to start the evolution of the heavy ion collision. This is motivated by fits to variables 
describing the collective behaviour of the fireball, such as Hanbury Brown-Twiss (HBT) radii and 
collective flow parameters. Pre-equilibrium flow is universal when describing a boost-invariant 
system with a traceless stress-energy tensor j9]. Since our model utilizes a 2+1 dimensional 
hydrodynamic evolution (see below), we are well-justified in employing pre-equilibrium flow in 
our description of heavy ion collisions. 
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Figure 2. Equation of state functions from HotQCD (solid green curve), s95p-vl (solid black 
curve), VH2 lattice-inspired EoS (pink curve), and the MIT Bag model EoS with a first order 
transition. 



Using the set up described above, we run VH2 2D+1 viscous hydrodynamics developed by 
Luzum and Romatschke [8j on the energy-density distribution that resulted from the above 
assumptions. At the end of the hydrodynamic evolution, we employ Cooper-Frye freezeout [10] 
with viscous corrections, and then start the UrQMD hadronic cascade [llj to obtain final particle 
distributions. 

After ensuring that the resultant yields and mean {{pt)) of the produced particles 
match those of experimental data, we employ our analytical framework to draw numerical 
comparisons between our model results and those from experiment. This is done by describing 
our modelling data using Chebyshev polynomials, and then calculating the /ndi between 
these polynomials and experimental data. The is calculated taking into account point-to- 
point correlated, uncorrelated, and scaling uncertainties (A, B, and C-type errors, as described 
in [12j). Performing this fit and calculation for a grid of initial state parameters, such as viscosity 
over entropy-density {rj/ s) and the maximum temperature of the evolving medium (Tcentrai)^ we 
are then able to do a global comparison of minima for a variety of soft observables: particle 
spectra, elliptic flow coefficients, and femtoscopic measurements. 

3. Data Analysis 

As mentioned in the previous section, randomized Glauber initial conditions for Au+Au 
collisions at ^/sNN — 200 GeV with an impact parameter (6) of 4.4 fm, shown in panels b 
and e of Figure [2} are used as an input for this study. The main goal of the analysis is two- fold: 
first, to study the differences due to inclusion or omission of pre-equilibrium flow. Second, to 
study the impact of using Ncoii vs. Npart scaling on the accuracy of describing experimental data. 
Four sets of modelling data were generated: two sets with pre-equilibrium flow turned on (one 
scaled with Ncoii, the other with Npart), and two sets with pre-equilibrium flow turned off (also 
one with Ncoii scaling and the other with Npart scaling). We varied Tcentrai (highest temperature 
of the medium) from 320 to 360 MeV, and ry/s from 0.001 to 0.36. The VH2 EoS was replaced 



i 



10' 



10 



' I ' ' ' ' I ' ' ' ' I ' ' " I " ' ' I ' ' ' ' I ' ' ' ' I ' 



_ ^ 4_— i>— I— i — 



• PHENIX data 

I I . . 1 I I ■ I ■ . I , . . ■ I ■ . . , I , 



50 100 150 200 250 300 350 400 

'^part 



V 



0.8 



0.4 



0.2 



' I " ' ' I ' ' ' ' I ' ' ' ' I ' 




• PHENIX data b 

I I . . , I I ■ I ■ . I , . . ■ I ■ . . , I , I I ■ I I . . . I . ■ . 



50 100 150 200 250 300 350 400 

'^part 



Figure 3. tt^, K+, and p+ dN/dy (panel a) and (px) (panel b) distributions from ^/snn — 200 
GeV 10-20% central Au+Au collisions, modelled using CHIMERA (collored filled symbols), VH2 
pre-errata data with Ts^itch^lSO (color curves), and experimental data from PHENIX (black 
and open symbols). 



by HotQCD HISQFIX EoS Cooper- Fry freeze-out was started when the temperature of the 
medium reached 165 MeV (Tfreeze-out), roughly corresponding to the temperature of chemical 
freeze-out derived from the thermal model [13j. UrQMD was started at Ts^itch = 140 MeV. 

To check for consistency, we ran a sample of Au+Au 200 GeV data without any centrality 
selections. We calculate yields per unit of rapidity (dN/dy) and (pt) as a function of Npart 
for the three most common particle species: protons, kaons and pions, and compare these with 
experimental data, as shown in Figure [3j Since the modelled and experimental data are in a 
reasonable general agreement, we can now start mapping out the precise correspondence between 
the experiment and modelling data with variable initial conditions. 

For each particle species, and for each of the three observables (HBT radii, spectra, and the 
elliptic flow coefficient), we obtain a set of data produced with parameters described above. An 
example of particle spectra, plotted together with experimental data, is shown in Figure [Ij The 
figure illustrates a typical model-data comparison done to show the validity of hydrodynamical 
codes, and demonstrates a fraction of the available parameters generated by CHIMERA. Note 
that each set of modelling data is fit using a Chebyshev polynomial (depending on the observable, 
the polynomial can have four or five parameters). 

For each experimental-modelling data pair, we calculate x^i quantify how well a given 
modelling curve describes a given set of experimental data points, as described in the previous 
section. Keeping all but two of the initial condition parameters constant, we are able to 
map out the x^/iidf space for particle spectra, flow, and HBT radii, and then use the most 
probable overlap region to determine the most likely initial states parameter values. In the 
study performed for this article, the x^/^df was mapped out for tt^ mesons, with experimental 
data obtained from both PHENIX [llj[l5] and STAR [16] [T7] collaborations. An example, shown 
in Figure [sj shows a calculation of one map point for each x^/ndf map (spectra, flow, and each 
of the HBT radii). Figure [5] calculations were performed with Npart scaled modelling data, with 
pre-equilibrium flow turned on, at 77/s=0.08, and Tcentrai=320 MeV. By eye, these parameters 
represent the data fairly well, but exactly how well can be judged by the x^/ndf shown in each 
panel of the figure. 




Figure 4. tt^ (panel a), (panel 6), andp (panel c) spectra distributions from CHIMERA and 
^/sNN — 200 GeV 20% central Au+Au RHIC data. 77/s is kept constant, while Tcentrai is varied 
from 320 to 360 MeV. Bottom sections of each panel show fractional agreement between data 
and a given set of parameters. Model data are produced with Ncoii scaling and pre-equilibrium 
flow off. 



Table 1. xV^df minima for each of the three experimental observables: tt^ meson spectra, 
flow (v2 coefficient), and the combined HBT radii. The minima have been computed with 
pre-equilibrium flow turned on and off, and with initial state scaling as the number of binary 
collisions (Ncoii) or number of participants (Npart)- 
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4. Results 

The results of the study are summarized in Table [l] and Figure |6] Table [l] shows the minima from 
X^/ndf distributions for all measurements currently available within CHIMERA made with tt^ 
mesons. The x^/ndf minima range from 4.25 for x^/ndf calculations of the combined HBT radii 
with Ncoii scaling and pre-equilibrium flow included, to 26.43 for the x^/ndf distributions derived 
from calculations for V2 without pre-flow and with Npart scaling. Trends, evident from the table, 
are as follows: x^/^df is significantly improved for V2 and HBT measurements by turning on 
pre-equilibrium flow and by using Ncoii scaling, while the reverse is true for pion spectra: the 
best x^/ndf value is for Npart scaling in simulations with pre-equilibrium flow turned off. 

Figure [g] shows x^/iidf space distributions for spectra, V2 and HBT radii of tt^ produced 
with Ncoii scaling (panels a, 6, and c) and Npart scaling (panels d, e, /). The distributions are 
drawn to display the most likely 77/s and Tcentrai values that would produce the experimental 
measurements. A global 77/s and Tcentrai value would be at the intersection of the x^/ndf minima 
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Figure 5. An example of x^/ndf calculations for tt^ mesons. Each panel shows a comparison 
between model and experiment for each of the available quantities: V2, spectra, and the three 
HBT radii (right side of the figure). Pre-equilibrium flow is turned off, and the initial distribution 
scales as Npart- Experimental data is shown in purple, and modelling data points - in salmon 
color filled circles. 



for all three experimental quantities modelled under the same set of conditions. However, this 
is not yet the case. In Ncoii-scaled data, 77/s and Tcentrai tend to higher values for spectra and to 
some extent for the HBT combined measurement as well, whereas for V2 the most likely 77/s is 
much lower 0.08). We also do not see a common minimum in the Npart-scaled data: there the 
best 77/s value for 7r+ spectrum is 0.32, 0.08 for V2, and 0.0001 for HBT. The best Tcentrai is the 
same for spectrum and HBT (320 MeV), but higher for V2. This indicates that the modelling 
part of CHIMERA is not yet completely physical. 

5. Summary and outlook 

In summary, we presented an outline and the first tests of an evaluative framework for comparing 
hydrodynamical modelling attempts of the Quark Gluon Plasma to the experimental data 
available from the existing heavy ion experiments. VH2 viscous hydrodynamical code was used 
as a base to generate modelling data used in this evaluation. The data comparison was done 
for 10-20% central ^/sNN — 200 GeV Au+Au collisions at RHIC, and aimed to study the effect 
of utilizing the initial scaling (Ncoii vs Npart)- In addition, we analysed possible effects of pre- 
equilibrium flow on the goodness of model's description of the data. The result of the comparison 
finds that Ncoii scaling might not be applicable when attempting to describe the spectra of the 
data derived from collision's bulk, but does a slightly better job than Npart scaling at describing 
flow and the HBT variables. We also found that pre-equilibrium flow significantly improves 
the goodness of fit for the HBT radii and V2, as expected, but also undermined the model's 
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Figure 6. Comparison of x^/ndf distributions for Ncoii vs. Npart scaling. The top row has been 
produced with Ncoii scahng, the bottom row scales with Npart- Panels a and d show x^/ndf for 
7r+ spectra, b and e - tt^ flow (V2), and panels c and / show the distributions for the combined 
HBT radii. Pre-equilibrium flow was turned on for the four panels on the right, and turned off 
for the left-most panels. 



description of spectra. 

However, these are initial results, generated primarily as proof of principle. Currently, we 
plan on implementing a more realistic initial-state energy-density distribution, vary the EoS 
used for the analysis, and move on to three-dimensional hydro-dynamical modelling. We are also 
performing Tgtart studies and studies to map out 77/s dependence on temperature. In addition, 
we are expanding our framework to analyse other collision centralities in Au+Au, other collision 
species and collision energies. The onslaught of new data from RHIC and the LHC provides us 
with a unique and exciting opportunity for a comprehensive and detailed study of the matter 
produced in these collisions. CHIMERA is an important tool for the Heavy Ion community, and 
we hope it would be used accordingly. 
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